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Abstract 

Using the operator representation of the Dirac Coulomb Green 
function the analytical method in perturbation theory is employed in 
obtaining solutions of the Dirac equation for a hydrogen-like atom in a 
time-dependent electric field. The relativistic dynamical polarizability 
of hydrogen- like atoms is calculated and analysed. 
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1 Introduction 



In stationary perturbation theory as well as in the time-dependent one, the 
method of using the Coulomb Green function has found wide application in 
obtaining analytical solutions of the Schrodinger or Dirac equations. The 
main advantage of this method is the possibility of obtaining the final results 
in closed analytical form or reducing results to the summation of rapidly 
convergent series. Therefore, by using the above-mentioned method one can 
avoid a lot of complex numerical integrations (see, for example, Makhanek 
and Korol'kov 1982, Zapryagaev et al 1985 and references cited therein). On 
the basis of the application of the connection between the problem of the 
four- dimensional isotropic harmonic oscillator and that of a hydrogen-like 
atom in eletromagnetic fields (see Kustaanheimo and Stiefel 1965), it has 
been proposed to establish a new representation of the Coulomb Green func- 
tion in the form of a product of annihilation and creation operators (this 
is called by us an operator representation). Operator representation of the 
Coulomb Green function is very efficient in application for the non-relativistic 
case (Le Van Hoang et al 1989) as well as for the relativistic one (Le Anh 
Thu et al 1994). The main elements of the algebraic method used (with 
the aid of the operator representation of the Coulomb Green functions) in 
Le Van Hoang et al (1989) and in Le Anh Thu et al (1994) are as follows. 
By using the above-mentioned connection, all operators of the algebra of 
the dynamical symmetry group 5*0(4, 2) can be found in the quadratic form 
of the annihilation and creation operators (see, for example, Kleinert 1968, 
Komarov and Romanova 1982). Therefore, the calculation method, based 
on the use of the algebra of 5(9(4,2), leads only to the use of the simple 
communtaion relations between the latter operators. The use of this method 
together with the operator representation of the Coulomb Green function 
essentially reduces the calculation process and provides for reducing rather 
complicated calculations of matrix elements with the Coulomb wavefunc- 
tions to a purely algebraic procedure of transforming the product of the 
annihilation and creation operators to the normal form. The advantage of 
the proposed algebraic method are found not only in the simplicity of the 
calculation process but also in the possibility of obtaining the final results 
in the summation of rapidly convergent series. In fact, some clever results 
are obtained in Le Anh Thu et al (1994) for the problem of calculation of 
the relativistic polarizability of hydrogen-like atoms. In this present paper, 
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we consider the problem of calculating the relativistic dynamical polarizabil- 
ity of the ground state of hydrogen-like atoms on the basis of application 
of the operator representation of the Dirac Coulomb Green function estab- 
lished in Le Anh Thu et al 1994. These calculations, besides their purely 
theoretical significance, are of great practical interest connected with recent 
developments in experimental investigations of multiply charged ions (see, 
for example, Zapryagaev et al 1985, Paratzacos and Mork 1979). However, 
the majority of accurate calculations has been done only for the static po- 
larizability of relativistic hydrogen-like atoms (see also Drake and Goldman 
1981, 1988, Johnson et al 1988). In our calculations the radiation corrections 
are neglected, taking into account the fact that this effect is small compared 
with the external field effect. Our results are directly generalized from the 
non-relativistic calculations (Zapryagaev et al 1985) and coincide with the 
results in the static limit (Le Anh Thu et al 1994). 

2 Equation in two-dimensional complex space 

The Dirac equation for a hydrogen-like atom in the field of linearly polarized 
light can be written as follows (h = m = c = 1): 



d 

-ici\ r-— 



+ (3r-Ze 2 + d - erx 3 (e ivt + e~ lut )j *(r, t) 



= ^ (!) 

where a\(\ = 1, 2, 3) and (3 are the Dirac matrices ; 9 and v are the amplitude 
and frequency of the external electric field respectively. Further on, we use 
the usual representation 

-(£)■ -Co)' Ho- : 

where and ^2 are two-component spinors and o\ (A = 1, 2, 3) are the 
Pauli matrices. 

The formal changes (see Le Anh Thu et al 1994) 
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reduce equation (1) to an equation describing the interaction between a "par- 
ticle" with complex coordinates £ s (s = 1, 2) and the external electric field. 
Here, in (2) summation is indicated by means of repeated indices. The scalar 
product of wavefunctions in £-space is defined by the relation 



J —oo J —oo J —oo J —oo 

(3) 

where £' s = ReC, s , = Im^ s . All operators appearing in (1) are henceforth 
considered to conform with the formal changes (2). Thus, the operator in 
the left-hand side of equation (1) is self-adjoint with respect to the scalar 
product of wavefunctions defined by (3). 

We will solve equation (1) by using the method of perturbation theory 
assuming the external electric field to be small. Its solution can be found in 
the form 

tf(r,i) = ¥°\r,t) + e¥ 1) (r,t) 

= ¥°\r) e- i£ot + 9u(r) e"**^-") + v (r) e - lt{£0+v \ (4) 

where \&(°)(r) is a wavefunction in the zero-order approximation, i.e. a solu- 
tion of the Dirac equation 



d 

—io>\ r- 



+ Pr-E r)¥ \r) = Ze 2 ¥°\r); (5) 
\ dx x J 

Eq is the energy in the zero-order approximation. By substituting (4) into 
(1) and taking into account (5), we then obtain the equations 



d 
dxx 



+ (3r- Ze 2 ^j u(r) - r(e - u) u(r) = ~ erx 3 ¥ 0) (r) (6) 



r^— +fir- ZeA v(r) - r(e + v) v(r) = -- erx 3 ¥°\r). (7) 
\ dx x J 2 
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Noting that equations (6) and (7) have the same structure we consider only 
equation (6) for the function u(r); then by replacing —v by v in u(r) we find 
the function v(r). 

Let us now present u(r) and in the form 



u 



(0) 



( I \ W ' * (0) = L / \ (0)1- 



(8) 



The substitution of (8) into (6) leads to the set of equations for ip^ and (p 2 



,(!) 



-i (^~^~ + X ) V ? 2 1) -^V ; '2 1) + r ( 1 -£ : o + ^)v ? i 1) -^ e Vi 1) = -^erx 3 </?S 0) (9) 

-i ^ A ^- + ^ ) +ik^ ) -r(l+e Q -v) V { 2 1) -Ze 2 ^ ) = ~erx^\ (10) 

where k = 1 + oa^a ; 1 is the orbital momentum operator. 
By using the transformations 

v4 0) = ~VT+T (^ - <?<°>) , J* = ±VT=T (jM + S (0) ) (11.1) 



and 



= -\y/l + e,-v(j™ _ gW) , ^) = Ivi-eo + f + 

(11.2) 

we find 



9 , Ze 2 . . 
x A - hl + wr (e - v) 



J™ +(k + —) gv 



(12) 



8 , Ze 2 . 
x x - hl-wr H (e - i/) 

OX\ Ul 



gw + U-^f) 



5 



= -rx 3 (A+FW + A-G<®) 



(13) 



where 



= yJl-( £o -u¥ 



A± = 



e 



^(l-e ){l-e + u) ± v /(i + eo )(i + e - I /) 



(14) 



After expanding the functions T and Q in power series of the eigenf unctions 
of operators L 2 and k, the operator k appearing in (12) and (13) becomes 
a c-number and, therefore, in order to solve equations (12) and (13) we can 
employ the Green function method established in Le Anh Thu et al (1994). 
In the next section we will show an example by calculating the dynamical 
polarizability of hydrogen-like atoms in the ground state. 



3 Relativistic dynamical polarizability in the 
ground state of hydrogen-like atoms 

Solution of equation (5) can be easily obtained by purely algebraic calcu- 
lations. In particular, we can find for the ground state the solution (see 
Komarov and Romanova 1985): 

V^o) 

Q {0) = 0, (15) 



where ujq = Ze 2 ; Eq = y 1 — (Ze 2 ) 2 ; are eigenvectors of operator a 3 and 
|0(wo)) is the vacuum state defined by the equations 

a 8 (u o )\0(u o )) = b a (u o )\0(u o )) =0, s = 1,2. 

Here the operators a s (u), b s (uj) are defined as follows 



cud£ s ) ' 



2 \, udtsj ' s v ' V 2 \, wd£ 
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where uj is a positive parameter (see Le Anh Thu et al 1994). The perturba- 
tion term in equations (12) and (13) thus has the form 

±A T rx 3 ^ = ±A T r £ «(Z 1 ^ 1 + v^Z 1)2 )|0(w )). (16) 
Here, we use the notation 

and Zi K are eigenvectors of (i) the square orbital momentum operator and 
(ii) the operator k. The structure of perturbation term (16) prompts us to 
find the solution in the form 



^ {1) = E ^\ s (1) = E <#> 

K=-l,2 



K=-l,2 

where and C^ 1 ) satisfy the equations 



5 , N 

X\- V 1 + LOT (£ - V) 



' dx 



A 



UJ 



_d_ 

dXy 



+ 1 - uir + 



-- A_^/\K\r et >Z lK \0(uj )}, 
Ze 2 



UJ 



-(e - f) 



UJ 



= -A + J\K\r £o Z lK \0{uj o )). 



From (19) it follows that 



(17) 



(18) 



(19) 



K — 



Ze 2 



r £ Z lK \0^o)) 



K — 



Ze 2 



9 , Ze \ 

x\ — h 1 -uor-i (e - v) 



'dx 



A 



UJ 



Substituting (20) into (18) we obtain 



8 Ze\ 
x\- Er - + 1 + ujt — - (e Q - v) 



(20) 



"dXy 



UJ 



d Ze 2 . 
xxtt~ + i-our-\ (e Q - v) 



' dxy 



UJ 
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+K*- 



Ze 



LO 



2x2' 



a) 



A. 



Ze 2 



K\ \ K — 



UJ 



8 Ze\ 
x\- Er - + 1 + ur — - (e Q - v) 



"dx> 



UJ 



r £o Z lK |0(^ o )> 



r £l % K |0(cu )>. (21) 



As mentioned above (see (2)), all operators in (20), (21) (though formally 
written for brevity via the usual coordinates x\ (A = 1, 2, 3)) are understood 
in the sence of the formal changes (2). These changes are equivalent to the 
transformation from r-space to £-space (Kustaanheimo-Stiefel transforma- 
tion, see Kustaanheimo-Stiefel (1965)): 

xx = SM*£t M = i,2) 

X = axg(fi) 

So we can rewrite the operators used in (20), (21) via the annihilation and 
creation operators as follows 



x\ 



d M(u) M + {uo) 



dx x 



1, 



r = —[M(u) + M + {u) + N(u) + 2] 
2cu 



where the notation 

N = afa s + bfb s 



M = a,6. 



M+ = atbt. 



Therefore, equation (21) has a more convenient form 



N 

M+- + 1 



Ze 2 



UJ 



(eo - v) 



N 

M+ + - + 1 
2 



Ze 2 



UJ 



(eo - v) 



-K 2 - 



Ze 2 \ 2 



to 



a. 



Ze 2 



K \ \ K — 



UJ 



r £ »Z lK \0(uj )) 



\, N , Ze2 , 
M+ — + 1 (eo-iy) 

z UJ 



r*°Z lK \0(uo))- 



(22) 



Here and henceforth, we omit for brevity the parameter uj in the expressions 
of the operators. Equation (22) has the same structure as the equations ap- 
pearing in the case of calculation of the static polarizability of hydrogen-like 
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atoms (see Le Ami Thu et al 1994). The only difference is in the perturbation 
term on the left-hand side of these equations. Therefore, we can by analogy 
find the solutions of equation (22), using the Green function operator which, 
according to Le Anh Thu et al (1994), can be established as follows: 



G lK = B lK , B lK = Y d s (N/2)M s , (23) 

lK 2 + N + M + M+' K 1 ' 



where 



d ° (n) = - n + 7 + 2-ZA/2 - (24) 

d in) = (-1)«(2 7 + 1) (" + ' + 1 > ! r(n + S+ l- 7 -ZeV2) 

x r(n + 2 + 7-ZA/2) 
r(n + s + 3 + 7-Ze%/2)' ' ' 

7 = -I - 1 + V/t 2 - Z 2 e 4 , e = e -v. 

Here, we note that the Green function operator (23) acts on the basis of 
states with frequency uo. However, the wavefunctions in the zero-order ap- 
proximation (15) have the frequency uj = Ze 2 . Therefore, in order to use 
the algebraic calculation we have to transform the wavefunctions (15) from 
the frequency u to uj using the unitary transformation (see Komarov and 
Romanova 1982): 

|$(w )) = U(u Q ,u)\&(u)), 

where 

U(u , u) = exp 1 1 In ^ [M(u) - M + (u)] } . 
This transformation can be reduced to the normal form 



fr( Wo , w ) = r ^exp(^M>))exp(i VM l n ^) 
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x cxp 



UJ — UJ 

UJ + uJq 



M(u) 



(25) 



Finally, by analogy with that was done in Le Anh Thu et al (1994) for the 
same equation, we find the solution 



<7 f 2 

\Q$) = A_J\k\U )r^G lK r e ^U(u ,u)Z lK \OM) + 

v uj ' 



+A + \\K\r 1 Gi K r~ 



M + — + 1 {eq-v) 

Z uj 



r £0 U(uj ,uj)Z lK \0{uj )). (26) 



The wavefunction \J-^) can be obtained by substituting (26) into (20). 

Let us now calculate the dynamical polarizability of the ground state of 
hydrogen-like atoms, the formula of which in £-space has the form 



a(u) 



This formula, after considering (8), (11), (15) and (17), can be written as 
follows 





erxs 




(l£(0) 


r|*(°)) 



K=-l,2 



where 



4u) u)J\k\ ... 



(27) 



By substituting the found solutions (20), (26) into (27) we find the expression 
for the positive frequency term of polarizability 



E I ti 



A\5 2 
B 



A 2 _B + 2A+A_S + 



TTK 



4u;(a> — uj ) 
(uj + UJ ) 2 



A+A_ + 



A\8 
B 



Auj 2 (uj - oj ) 2 A\ 
(uj + uj )* ~B 



(28) 
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Here, we use the notations 

Ze 2 



B = K 

UJ 



UJ + UJ UJ 



HZ m = (0(^)|Z lK f> + (cu o ,cu)MV £o ^G K r £o ^(M + ) m f>(cu o ,cu)Z lK |0(cu)). 

A similar expression for the negative frequency term can be obtained by 
replacing v by —v in formulae (14), (28), (29). 
By using the algebra of operators M + , M, N: 



M, M + \ = N + 2, [M,N + 2] = 2M, [N + 2, M + \ = 2M + 
and the relations 

N(M + ) n Z lK \0) = 2(n + l)(M + ) n Z lK \0), 
M(M + ) n Z lK \0) = n(n + 21 + l)(M + ) n ^Z lK \0), 

" ZlM ~ (2ujyT(-p) h s\ ( S + 2/ + l)! (M } ^ |0) ' 
we algebraically obtained the explicit form of the term H£ as follows: 



r(g + go -7)r( g + 3 + g o-7) ^ r d s (p-s + l) p ° u t(p-s" 
r(g + 5 - 7 -m)(g + 3)! ) (P " *)< ST U * , 



x 



w + cjq/ r(t+ 1 +e + 7 -n) (t + 3)! ' 1 ' 
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Direct calculations show that all power series appearing in the term H£ m 
are rapidly convergent. The high convergency of these power series is di- 
rectly related to the expansion (23) of the Dirac Coulomb Green function 
which, in fact, is established on the basis of harmonic oscillator wavefunc- 
tions. Moreover, the expression (28) for the polarizability contains only H* m 
with n, m — 0, 1, the calculation of which needs only some first terms in the 
summation over p. For example, for frequency v less than the fourth reso- 
nance frequency (relative to the transition from ground state to 3P3/2 state) 
the contribution of the terms with p — 0, 1 in H* m (n,m = 0, 1) is about 
98 %-99 % for all Z < 137. 

In figures 1(a) and 1(b) the dependence of relativistic polarizability on the 
external field frequency is given for hydrogen-like atoms with Z = 50. The 
dotted lines correspond to the non-relativistic limit case. Figure 2 give the 
same for Z = 100. 

Let us now consider the non-relativistic limit, i.e. take into account only 
the first term in the expansion in the power series of Ze 2 . For this case, we 
can effect summation over the variables t, q, s in the expression of H£ m and 
then have H£ m in the form of hypergeometric functions. Consequently, we 
have 



2 n /i 3 e 2 J 1 / 1 1 ffi-V 

> ±v) ~ (z^ + iri(2-i)( 3 -i) 2i H 5, / ^ U+r. 

5u 2 / 1,1(11-1 

+ 7 , ,1 1T2F1 6,3 --,4-- ' 



Gu + l) 2 (3-i) ^ ' /i' /x'V/x + 1 



where [i — yl ± 2u/(Ze 2 ) 2 . This expression coincides with the well known 
result obtained by Vetchinkin and Khristenko (1968). 

By putting v — into (28), (30) and (31) we thus find the formula for the 
relativistic polarizability. It is easy to see that for uj = uj the formula (28) 
contains only the term Hq . Taking into account the formula (Prudnikov et 
al 1981) 

we can lead Hq to the form 
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T(e - 7) r(£ -7 + 3) r(e + 7 + 1) r(e + 7 + 4) 



x 



E 



E 



d s ip - s + 1) 



(32) 



P =o r ( £ o - 7 - p) s = (P - s ) ! (P " s + 3 ) ! r ( £ o + 7 + 1 - V + s) ' 

By substituting (32) into (28) we obtain for v = the relativistic static 
polarizability 



e 2 (e + l)(2e + l)(4eg + 13e + 12) e 2 (e - 2) 2 r(e + 7 + 4) r(e - 7 + 3) 



36c^o 



36^ 4 r(2e )r(-eo - 7)r(l - e + 7) 



x 



~ r(fc- g0 -7)r(A:-e + 7 + l) , ^ T(q - e + 7 + 1) 



Lfc=0 



fc! (A; + 3)!(A;-e + 7 + 3) 



q=l 



(9 + 3)! 



^ r( g - £ - 7 + 2) ^ r(s - £ - 7) r( s - £ + 7 + 3) 

X r(g-5 + 7 + 4)^ ' s !r((s- £o -7 + 3) 



where 7 = — 2 + y4 — (Ze 2 ) 2 . This result absolutely coincides with the result 
obtained in Le Anh Thu et al (1994) (see also Barut and Nagel 1976). 



4 Conclusion 

In conclusion we would like to note that the magnetic field effects, as a rule, 
should be taken into account for a detailed investigation of the behavior of 
a relativistic atom in the field of linearly polarized light. These effects can 
be neglected only in the non-relativistic limit. The above method proposed 
with the use of the operator representation of the Coulomb Green function 
can also be employed, for example, in calculating magnetic polarizability. 
Consideration of magnetic field effects leads only to an enormous number 
of calculations, which are more complicated in comparison with the above 
calculations but could be done analogously by analytical methods. Moreover, 
we hope that our algebraic method would be helpful when considering the 
problem of an atom in a quantum field, in particular in calculating the Lamb 
shift of multiply charged ions, which is of great interest and has been widely 
investigated recently (see, for example, Snyderman 1991). 
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